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Abstract 

We demonstrate that the corrections to the classical Kelvin image theory due 
to finite electron screening length A, recently discussed by Roulet and Saint 
Jean, Am. J. Phys. 68(4) 319, is amenable to an exact closed form solution 
in terms of an integral involving Bessel functions. An improper choice of 
boundary conditions is rectified as well, enabling also a complete solution for 
all potentials - both inside and outside the metal surface. 

I. INTRODUCTION 

The theory of image charges over a perfect metal surface has an early history that dates 
back to Lord Kelvin in 1848 ||l|. It is a subject covered in most textbooks on electrostatics 
0, 0, as well as an important topic in modern research Hence it is of considerable 

importance in the undergraduate curriculum. In a recent article in this journal, Roulet 
and Saint Jeans (RSJ) discussed the corrections beyond the classical Kelvin theory by 
considering the effect due to a finite electron screening length A in the metal. This topic 
has indeed hardly been addressed in any textbooks, except for more advanced texts such as 
Mahan , but only in the more general theory of electron screening and response functions. 
A simplified consideration based on chemical equilibrium such as employed by RSJ P] should 
deepen the understanding of classical Kelvin image theory as well as stimulate advanced 



students to consider the more comprehensive treatment based on many-body theory and 
hnear response Unfortunately the otherwise excellent exposition of RSJ [|| is plagued 
by an improper choice of boundary conditions, and by their inability to obtain a complete 
solution for the potential in the metal dictated by a Helmholtz type partial differential 
equation (PDE), which as we shall see is in fact separable, see Appendix I and II. They 
resorted to a perturbative solution that is not valid for large screening lengths - which is in 
fact a trivial limit. In this paper we shall rectify this and consequently provide the complete 
solution for the potentials both inside and outside the metal surface. Naturally our results, 
being exact will recover the expected classical theory both in the case of small screening 
length (A 0), as treated by RSJ 0, and the opposite case of large screening length 
(A oo), when the metal becomes ineffective and the classical Laplace potential ensues. 

This paper is divided into four sections. In section I we shall review the classical Kelvin 
image theory, solved by the separation of the Laplace equation in cylindrical co-ordinates. 
This detail solution is nowadays not commonly taught in the undergraduate curriculum 
1^. In section II we shall demonstrate that by using the same technique, the Helmholtz 
equation (Eq.(22) of ref for the screened potential inside the metal can also be solved. 
More importantly we shall discuss the proper boundary conditions for this problem. The 
standard conditions follow from the Maxwell equations 0: {div D = p and curl E = 0) 
on the surface edge {z = 0). However RSJ imposed an artificial condition based on a 
strict compliance of the final surface charge density (here denoted as cr(r)) with the limiting 
classical surface charge density (here denoted as cTo(r)). Here we shall argue that the limiting 
surface charge density o"o(r) should only be a consequence of the complete theory and not 
a precondition. As a result a{r) should be A dependent, while cro(r) is manifestly not. 
The RSJ choice of boundary condition, as it will be shown is also amenable to a closed form 
solution (see Appendix II), hence their perturbative approach is unnecessary. It is interesting 
to note (see Appendix II) that the exact solution to their problem in fact diverges as A^ 
for large A (see Appendix II) and thus the perturbation method is inoperable in principle, 
as it can never converge to the exact solution. We believe this essentially non-perturbative 



feature is an artifact of their boundary condition which is unhke the solutions presented here. 
Moreover, although the departures between our theory and theirs only become apparent at 
O(A^) for small A, their theory further leaves the outside potentials undetermined and will 
require additional assumptions to complete. We shall show that our solution to the problem 
is complete and can be expressed in terms of integrals involving Bessel functions of the 



type similar to and much studied in electromagnetic propagation 0, |[T0|, beginning with 



Sommerfeld's classic work of 1909 |]TT|. In section III shall present the exact solutions and 
in section IV we shall discuss the effects on and hence corrections to the classical image 
potentials outside the metal surface, which were omitted by RSJ 0. Here we shall also 
obtain the corrections to the surface charge density as discussed above. This correction 
integrates to zero charge as we shall see (Appendix III). In section V we shall replace the 
metal in our theory by a dielectric with e > cq. The reader may like to note that all the 
results of this paper are not new. They have in fact been derived before (unknown to the 
author) by Newns |12| and more recently by Krcmar et al (see the end notes). As a result 
of our investigations, the conclusions reached by RSJ for the case of large A//i>>lasina 
semiconductor is subtle as naively the limit A ^ cxo in fact approaches the classical results 
as required for a dielectric. A more careful comparison will require investigating the case 
X » h under the condition that 2;, r >> A which is indeed a non-trivial limit in our theory. 



More detailed numerical calculations can be found in ref. ||T3[. We shall conclude with a 
brief discussion of time-dependent effects, such as the case of an oscillating charge, by which 
ingenious experiments based on the classical skin effect could be used to test our predictions. 
These experiments could be designed and possibly be inspiring to an undergraduate class. 



II. CLASSICAL IMAGE THEORY 

In this section we shall adopt the more common convention in which the physical charge 
q is introduced above the metal surface a.t z = h > and where the rest of the metal is 
defined for 2; < 0, see Fig 1. This is opposite to the case of RSJ but is similar to most 



textbooks such as Jackson [0] and Landau et al |^. It is well known that Kelvin image 
theory provides the solution for the potential outside the metal as: 

0(r, z) = , ^ , ^ (2.1) 

A-ne^Jr'^ + {h - zY A-ne^Jr'^ + {h + zY 

The second term is the so called image potential due to the fictitious charge —q whose 
addition enables the boundary condition 0(r, 0) = to be satisfied, see Fig 1. This is a 
Dirichlet boundary value problem whose solution Eg (|2.1| ) is justified heuristically in most 
texts, see for example 0. It is perhaps useful in a later undergraduate course to actually 
verify that Eq ( p.l|) is indeed a solution of the Laplace equation. This can be demonstrated in 
the following way. The Laplace equation in cylindrical coordinates with azimuthal symmetry: 

has solutions that are given in terms of Bessel functions P , |14 



roo 

r,z)= / f{k)e~''''Jo{kr)dk for z > h, 
Jo 

roo 

/ g{k)sinh{kz)Jo{kr)dk for < z < h. (2.3) 
Jo 



The coefficients of these expansions namely f{k) and g{k) are to be determined by appro- 
priate boundary conditions as we shall see. Note that the second expansion has been chosen 
with the sink function in order to satisfy the required boundary condition on the metal 
surface z = 0. We have two coefficients and thus we need two more boundary conditions to 
determine the solutions. As is well known from most texts [0, 0, |1^, for electrostatics, 
these are derived from the two Maxwell equations: 

curl E = and div D = p, (2.4) 

which upon integrating around an infinitesimal surface volume imply the continuity of the 
potential, say at z = h: 

<P%{r,h) = 4>l{r,h), (2.5) 



and the discontinuity of the slope: 

eo ^ oz oz ^ z=h 

where the charge density p consists only of the true charge and does not include any induced 
charges. The delta function in Eq.( p.6|) is normalized for: 

POO 

2n / r6{r)dr = 1. (2.7) 
Jo 

Now the continuity condition Eq. ( p.5|) leads to: 

f{k) = g{k)e''''sinh{kh), (2.8) 

while the discontinuity condition Eq. ( p.6| ), with the use of the orthogonality property of 
the Bessel functions jr5| : 

POO I 

/ rJo(kr)Mk'r)dr = -6{k-k'), (2.9) 
JO k 

leads to: 

^ - g{k)cosh{kh) + f{k)e-''\ (2.10) 



The straightforward solution of the simultaneous Eq. (|2.8|) and Eq. (|2.10|) leads to: 

f{k) = sinhikh) and 
27reo 

9{k) = -^e-'\ (2.11) 

That these solutions lead to the classical Kelvin image potential Eq. (|2.1|) follows readily 
from another student exercise in the form of a mathematical identity for Bessel functions 

g: 

1 r'^ 

e-^\^Kjo{kr)dk. (2.12) 



y/r^~+~h^ Jo 

Notice that the final solution depends on the boundary conditions both aX z = h and at 
z = 0, and not just on the latter alone as the heuristic argument seems to suggest. 



Finally to conclude this section we shall derive the classical induced surface charge density 
cro(r). This is simply obtained from the required boundary condition |T^ : 



cro(r) = -eo ^ 



dz 



^ e'^\''\kJo{kr)dk = - J (2.13) 



2 = 



27reo7o ' ' 27reo (r2 + /i2)3/2 ^ 



a result easily derived using Eq. (|2.12|) . We remind the reader that this boundary condition 



is obtained from an infinitesimal small surface volume integral of the Maxwell equation: 

div E = -, (2.14) 
eo 

in the same way as Eq( p.6| ) where the total charge density pt contains all charges, true and 



induced [|I4[. In the section IV, as with RSJ 0, we shall derive this surface charge density 
by integrating the bulk interior charge density p°(r, z) as we shall see. 

III. BULK CHARGE DENSITY - AN EXACT SOLUTION 

The (Thomas-Fermi) screening theory modifications for the potential inside the metal 
(pin is given by the Helmholtz equation 0: 

^+-^(4)*"=% t« ^<°. (3-1) 

oz"^ ror^or-' 
where the screening length A is given by: 

This is expressed in terms of a thermodynamic derivative of the chemical potential with 



respect to the density. We note that this theory for A given by eqn(|3?T|) is in fact quite 
general, it should also be valid in the case of a dielectric or a liquid surface, upon replacing A 
by the appropriate Debye-Hiickel screening length. Significant modifications to the Thomas- 
Fermi theory will come about only when A becomes comparable to the Fermi wavelength 
\f which we shall briefly mention in the conclusion. In Appendix I, we shall show that the 



solution of the Helmholtz equation Eq. (|3.1| ) is also expressed in terms of Bessel functions. 



Hence the complete solution of boundary value problem specified by Eq. (|2.2| ) and Eq. ( pJ] 
is now given by: 



■'>V' 5 



poo 

r,z)= / f{k)e~''''Jo{kr)dk for z > h, 
Jo 

^^{r,z) = / {gi{k)e^' + g2{k)e~''')Jo{kr)dk for 0<z<h and, 

JO 

i-co U 

0,„(r,^)=/ g^[p^)e^^^ — J^{kr)dk for ^ < 0; (3.3) 

where (31 = k"^ + V-^^j see Appendix I. Note that the non- vanishing of the potential on the 
surface {z = 0) for general A requires that we construct solutions for 0< that contains two 
(in general) unequal coefficients gi{k) and g2{k). We have already discussed the boundary 
conditions aX z = h. The continuity condition Eq.(|2.5|) leads to: 



f{k)-g2{k)=g^{k)e''\ (3.4) 

while the discontinuity condition Eq.( p.6|) leads to: 

9i{k) = T^e-'=^ (3.5) 

Two more boundary conditions are required to obtain the solution. These are in fact ob- 
tained in the same way as Eq.( |2.5|) and Eq. (p.6|) , only that now we have continuity of the 
potentials and their slopes on the surface 2; = in the absence of true charges on the surface. 
Thus from: 

0in(r,O) =0<(r,O), (3.6) 



we obtain: 



9i{k)+g2{k) = ^gM, (3-7) 



and from: 



dz 



90< 



=0 dz 



(3.8) 

2=0 ^ ' 



we obtain: 



(3.9) 



The problem is now completely specified by equations: (|3.4|),(|375|),(p.7|) and ( p.9| ). A si- 
multaneous solution of these equations leads to a complete solution for all the potentials, 
hence: 



^>[r, z) 



27reo ^0 

q roo 



kcosh{kh) + j3ksinh{kh) 



27reo ^0 
q roo 



k + Pk 
kcosh{kz) + (3ksinh{kz) 



e ^'^jQ{kr)dk for z > 



k + Pk 



-kh 



Jo{kr)dk 



e-'=Vo(A;r)dA; for z < 0; 



for < z < h and, 

(3.10) 



2neo Jo k + Pk 

In the Appendix II, we shall present in the same way the exact solution for the RSJ problem 
which consist of Eq. (|3.1| ) only and the boundary condition that a{r) = cro(r), which however 
leaves the outside potentials unspecified. 



IV. SCREENING CORRECTIONS 

Naturally it would be meaningful to compare these potentials with the classical image 
potentials in Eq.( p.3| ) and Eq.( |2.1lD . To do this each coefficient in Eq.( |3.1Cl| ) is rationalized 



by multiplying with {k — Pk)/ {k — Pk) and after some simple algebra, we have: 

(t)< = <p< + 0^ and 0> = 0> + 0\ (4.1) 

where (fP^ and 0> are the classical potentials given in section II. The screening corrections 
are given in terms of the A dependent potential: 

'^'^'"'^^ = ^ ^~(\/^' + - k)e-^^''^''hUkr)dk. (4.2) 
We can easily derive that for small A: 



^^(r.z] 



r ,-Hn.^)kUkr)dk = (4.3) 



A-o 27reo^ ' ' 27reo [r^ + (/i + 2)2]3/2 ' 



For large A we use an expansion of the square root:wl + l/a;^ ^ 1 + l/(2a; ) to obtain: 



6^(r, z] 



(1 r ^-Kh+z) j^^j^^^^^ ^ _A ^ (4.4) 



A-^oo 47160^0 " -^v- /-- 47reo ^r'^ + (/i + z)2^ 

which cancels the image potential in Eg . (12.11 ) as required. Using the same manipulations we 
obtain the inside potential as: 



27reo 

from which the bulk charge density: 



^"(WP + - k)e-^''ey^'^^kJoikr)dk, (4.5) 



p{r, z) = -^0.„(r, z) = ^l + ^je" e V A;Jo(A;r)rffc (4.6) 

is obtained. In particular we have the limit: 

in agreement with RSJ 0. However it is easy to show that the higher order terms differ, 
see Appendix II. Moreover unlike RSJ, we have no difficulty with the limit of large A which 
again by an expansion of the square root leads to: 



r, z 



which is of course the correct Laplace potential for the case when the metal is ineffective, 

since the Helmholtz equation Eq (p.l|) now reduces to the Laplace equation. As promised in 
section I, we shall evaluate the surface charge density: 

a{r)= p{r,z)dz = ao{r) + a^{r), (4.9) 

J —OQ 

where 

a^(r) = ^ / , k^Jo(kr)dk, (4.10) 



'k^ + 1/A2 

which is the new term in our theory. This quantity integrates to zero charge, (see Appendix 
III), thus the total charge remains —q as before. Hence the form do is not the only (unique) 
surface charge density that integrates to a total charge of —q. We note that while the limit 



A ^ oo appears straightforward, this is deceptive. Attempts to calculate the next order i.e. 
0(1/A^) corrections lead to divergent integrals that require careful treatment. We shall not 
discuss this exercise here. Suffice to say the corrections to the classical potentials due to 0^ 
and (pin are to smear the image charge from a point charge to a charge distribution whose 
weight vanishes as A increases. This charge distribution can be derived from the results 



presented here and has been worked out by Newns [0 in terms of an order two Bessel 
function. 



V. SEMICONDUCTOR SURFACE 

The extension of our theory to a semiconductor surface, in which A is finite, is now rather 
straightforward. The only significant modification is that the boundary condition Eg. ( p.8|) 
is now replaced by: 



dz 



^=0 ° dz 



2=0' 



(5.1) 



where e > eo is the dielectric constant of the material. We need not repeat the details here 
and present merely the solution for 0j„ which now takes the form: 



[r,z] 



'^^ -kh 



e-''''Jo{kr)dk for z < 0, 



(5.2) 



27reo ^0 k + ePk 

where e = e/eo- The integral can be analyzed in both limits as usual. For small A we have: 

gA h 



,2/A 



27reeo [r2 + /i2]3/2 
which is analogous with the metallic case. For large A however: 
-^/;e-VVo(^r)..^^^ 



{z < 0), 



(5.3) 



+ ejeo Jr^ + {h - 



{z < 0), 
(5.4) 



which is the classical result. Recall that the potential inside a dielectric is equivalent to an 



image charge q given by ||T4| 



replacing the real charge q at the point h outside the surface. Hence for large X >> h, 
classical theory is recovered and as a matter of fact it is for A << /i, as shown in Eg. ( p. 3D 
that screening creates departures from classical theory. Thus for dielectrics it is the case of 
small and finite A that leads to corrections to the Kelvin theory. Finally we shall briefly 
mention the case of an oscillating charge. The results here being static will lead to a 
zero frequency contribution to the classical finite frequency skin effect. The propagation of 



electromagnetic waves in the substrate will now be dictated by a complex wave vector [17 



= iuj'/c' + 2t/5' - 1/X' (5.6) 



to which we have added the screening length A. At low frequencies, [|T^, the skin depth S 
which varies as the inverse square root of the frequency dominates. Hence by extrapolating 
skin depth measurements to zero frequency, the results presented in this paper may be 
detectable |]l8l . This could be a novel undergraduate experiment for physics and engineeering 



classes. 



VI. CONCLUSION 

In conclusion we have presented an exact solution for the corrections to the classical 
Kelvin image theory of electrostatics previously discussed by RSJ in this journal. Some 
inadequacies in their analyses are rectified and a complete solution for all the potentials 
are obtained in closed form. We found that some care needs to be exercised with regard to 
statements about non-classical behaviour which are only valid in the case of finite A, whereas 
for small and large A we have shown that the theory reduce to the standard textbook analysis. 
Nevertheless our theory shows that the no n- vanishing potentials 0^ and for any finite 
A > in the case when the charge is right on the metal surface i.e. for /i = as opposed to 
the classical theory, have an essential role for surface chemistry that should be noted in all 
textbooks ||19| . 



A. Notes added: 



After this work was completed we were grateful to be advised by Gabriel Barton and 
Bernard Roulet that the results of our paper, with the exception of the proofs in the ap- 
pendices, have in fact been derived some thirty years ago by Dennis Newns [1^ apart from 
minor differences in our definition of the induced potentials. Also unknown to the author, 
the results have also been rederived recently by Krcmar et al ||TB|. The results of this paper 
are in agreement with Newns and Krcmar et al The reader might be interested to know that 
the latter authors have also considered the case when the charge q becomes immersed inside 
the material, using extensions of the methods presented here. Neither Newns, Krcmar et al, 
nor the present author have considered the limit when \p becomes significant. The author 
is further indebted to Marshall Stoneham for pointing out that the solution in the limit of 
large >> \tf has been considered by John Willis only within a static charge ap- 



proximation. This identified the significance of Fermi surface effects (omitted in this work) 
which lead to Friedel type oscillations in the density near the surface. No doubt a more 
sophisticated theory such as one using a density functional approach will be needed to treat 
all the various issues. 

B. Appendix I - solution of the Helmholtz equation 



The Helmholtz equation Eq.f p.l]) in azimuth symmetric cylindrical coordinates can be 
easily separated |16| by the ansatz (f){r,z) = f{r)g{z), such that: 

which has solutions: f{z) = e^^, where fc^ is the separation constant. Note that we have 
2; < so that the other solution f{z) = e~^^ is exponentially increasing and cannot be 
admitted by the boundary condition for 2; — > — 00. The radial equation now takes the form: 



where: a| = (fc^ — ja). For > this is the differential equation for Bessel functions 
of order zero. For = the solution goes as logr which is inadmissible as is the case 
where < since Bessel functions of imaginary arguments namely /o('") and -ft'o(^) are 
also inadmissible. Thus the general solution of Eq.( p.l| ) is a superposition of solutions using 
an expansion coefficient g^lk) given by: 

POO 

(f)in{r,z)= g^{k)e^' Jo{akr)dk for z <Q. (6.3) 

By a straightforward change of variable x = we easily obtain the solution given by the 
last of Eqn.(U). 

C. Appendix II - exact solution of the RSJ problem 

The solution for the RSJ boundary value problem can be obtained using the results of 
Appendix I. The RSJ potential and hence the bulk charge density is given by: 

p{r, z) = -^ g^lf3,)e^>^^ j^lkr)dk for z < 0. (6.4) 
Jo /3k 

The RSJ boundary condition requires that the integral: 

/O £q roc 

dzp{r,z) = -— / gz{l3k)-^JQ{kr)dk 
-oo A JO Pk 

= ao{r) = -— r e-^^kJQ{kr)dk (6.5) 

271 Jo 



hence we obtain: 



9sWk) = ^f3le-'\ (6.6) 



Note that while this determines the inside potential in closed form: 



RSJr^^ 2) = ^ e~^^e^^'l3kkJo{kr)dk, (6.7) 



in 



27reo "'0 

it leaves the outside potentials undetermined since Eq. ( p.4| ) and Eq.( p.5|) now become two 
equations with three unknowns. Moreover Eq.( |6.7|) is now divergent as for A — oo and 



does not reduce to the classical Laplace potential as required. It is interesting to compare 
Eq.( |6.7| ) with Eq.( |4.2| ). The former is essentially non-perturbative whereas the latter may be 
obtainable perturbatively if some care is exercised in treating the asymptotic integrals, see 
the remarks at the end of section IV. It is unclear and it would be interesting to speculate if 
this feature of the RSJ solution has experimental significance such as when the total surface 
charge can be maintained constant, i.e. cro(r) = const. In this case, following the Eq. (|6.5| ), 
we can show that the potential now goes as Ae~^/^. Nevertheless it is interesting to compare 
the small A expansion for the two solutions. We note that the difference between Eq.( |6.7|) and 
Eq. ( [4.5|) is a negative term: 

(f)Ur, z) = ^ e-^V'=^(A - k)kJo{kr)dk. (6.8) 



27reo "'0 

For small A we change variable to k = x/h and carry out the expansions for the square root 
and exponential terms via: 



'l + A'^l + lA' + ... and (6.9) 



h ' 2' h 



ef (v/^^) ^ (1 + 1^ + ...). (6.10) 



Inserting these into the integrands and evaluating the integrals we have the limiting expres- 
sions to 0((^)^) as: 



in agreement with RSJ while our solution has the expansion: 

qXe'/^ h N 3A {2h'^ - 3r^) _ X {2h^ - r^) i 

^' 27reo (/i2 + r2)3/2 + + ^' (/,2 + ^2)2 - I (h^ + r^) " " "J ' ^^'^^^ 

which differs from the latter at 0((^)^) with an additional term . 

D. Appendix III - proof that the integral of a^{r) vanishes 



The proof is straightforward and it makes use of a simple trick involving a well known 
property of the Bessel functions. We need to obtain, denoting 1/A by a, the integral: 



Q{a) = q rdr dk =A;Vo(A:r) 
"'0 Jo y/k^ + 

foo foo p—kh 

lim 



roo roo g 

q rdr dk ^ fcVo(/3r) Jo(fcr), (6.13) 

./o Jo V + a 

which follows from the property that lirrix^o Jq{x) = 1- Making use of the orthogonality 
property Eg. ( p. 91) we now have: 

Q{a) = limp_,o q dk —6{(3-k) 
Jo V + P 

= lim^^o q- 



= for « > 

= q for a = (6.14) 

respectively. Although Qia) is a discontinuous function at a = 0, the interchange of the limit 
with the integrals, while not rigorous, may be justified by the convergence of the integrals 
in Eq.dO^) for all a > |T|. 
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FIGURES 




Fig. 1 Cylindrical coordinates for the charge q placed at a distance z=h outside 
the metal or dielectric surface at with its image charge -q at z— h. The 
.solution for the potentials at any point P(r.z) with screening corrections is given in 



